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Abstract 

We establish a one-to-one correspondence between (1) quantum 
teleportation schemes, (2) dense coding schemes, (3) orthonormal 
bases of maximally entangled vectors, (4) orthonormal bases of unitary 
operators with respect to the Hilbert-Schmidt scalar product, and (5) 
depolarizing operations, whose Kraus operators can be chosen to be 
unitary. The teleportation and dense coding schemes are assumed to 
be "tight" in the sense that all Hilbert spaces involved have the same 
finite dimension d, and the classical channel involved distinguishes d 2 
signals. A general construction procedure for orthonormal bases of 
unitaries, involving Latin Squares and complex Hadamard Matrices is 
also presented. 

1 Introduction 

Teleportation and dense coding are two processes, which stood at 
the beginning of modern Quantum Information Theory. They both 
demonstrated radically new features of quantum information as op- 
posed to classical information, in that both would be impossible with- 
out the assistance of entangled states. Indeed, the attempt of using 
the properties of a classically correlated system shared by sender and 
receiver to improve the transmission rate of a classical channel can 
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easily be seen to be hopeless. But this is precisely what happens in 
teleportation and dense coding, and dramatically so, because without 
entanglement assistance teleportation, i.e., the transmission of quan- 
tum information on a classical channel, would not only be less efficient, 
but virtually impossible. 

In the original papers [ BW , BB] the new possibilities were demon- 
strated by giving an explicit example, based on qubits. It was clear 
early on that extensions to systems with higher dimensional Hilbert 
spaces were possible, not only to powers of 2, by running the process 
several times, but to any dimension 2 < d < oo |BB| . 

The task set in this paper is to do this systematically, and to 
classify all schemes for teleportation and dense coding. There are 
several reasons for doing this. The first is, of course, to take these 
miracle machines apart and to analyze what makes them work: what 
is the mathematical structure one really needs to set up such a scheme? 
For the present author one motivation of this kind was to understand 
the surprising observation that each of the published teleportation 
schemes also works as a dense coding scheme, and conversely: sender 
Alice and receiver Bob merely have to swap the equipment they use. 
An attempt at a direct proof of this failed, and indeed, as discussed 
below, the statement fails in general, but is true in the special case of 
"tight" schemes. 

The second reason for attempting a complete classification of tele- 
portation schemes is more practical. In spite of amazing progress in 
recent years, experiments in quantum information processing are still 
quite difficult. Hence, for realizing a teleportation scheme it is use- 
ful to have a systematic overview of the options, before going on to 
find the one which is the easiest to implement. This also goes for 
approximate realizations. And in order to find feasible approximate 
teleportation schemes it is probably once again necessary to under- 
stand the manifold of exact realizations. 

The aim of determining all schemes is not quite achieved in this 
paper, in two respects. Firstly, we will only look at the case when 
dense coding and teleportation are realized optimally with minimal 
resources, in the sense of Hilbert space dimensions and number of 
distinguishable classical signals. As in the well-known qubit case, this 
means that an entangled state between systems of the same dimension 
d as the input systems is used, and the classical channel distinguishes 
d 2 signals. That is, the classical capacity of the quantum channel 
is exactly doubled by dense coding, and teleportation requires twice 
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as much classical channel capacity as the quantum capacity of the 
channel set up by this scheme. We will call schemes with these di- 
mension parameters tight. As mentioned above, for these dimensions 
the symmetry between teleportation and dense coding holds perfectly. 
Classifying all schemes beyond the tight case appears to be more diffi- 
cult because there is too much freedom, which cannot be parametrized 
in a simple way (see, however, pDfl ). 

The second respect in which this paper falls short of a complete 
classification is that we can only reduce it to another "standard" prob- 
lem, namely the construction of orthonormal bases of unitary opera- 
tors with respect to the scalar product (A, B) ^ d^ 1 tv(A*B). In the 
last section we provide a fairly general construction for such bases. 
However, even this construction has to rely on other well-known but 
not completely classified combinatorial designs, namely Latin squares, 
and complex Hadamard matrices. This suggests that a complete con- 
struction procedure for all unitary bases would be at least as difficult 
as a complete classification of Latin squares or Hadamard matrices, 
and hence hardly a promising task. 

The paper is organized as follows: in Section 2 the Main Theo- 
rem is stated: an equivalence in the tight case between teleportation 
schemes, dense coding schemes, orthonormal unitary bases, bases of 
maximally entangled vectors, and so-called unitary depolarizers. Ba- 
sic consequences of the Theorem are discussed. Section 3 contains 
the proof, divided into subsections, each devoted to some implication 
in the big equivalence. In writing the proof an attempt was made 
to include also simple steps explicitly, and to make as transparent as 
possible why the tightness condition is crucial. Finally, in Section 4 we 
present the "Shift and Multiply" construction of unitary bases, which 
are then classified in terms of Latin squares and Hadamard matrices. 



2 Main result 

In order to state our result we use the following notation and termi- 
nology: When H is a Hilbert space, we denote by B(Tt) the space of 
bounded linear operators on Tt. A channel converting quantum sys- 
tems with Hilbert space 7ii n into systems with Hilbert space TC out is 
a linear operator T : B(7i ou t) — > B(7i m ), which is completely positive 



[Da, Pa | and normalized as T(l) = 1. A (discrete) observable F on 



7i over an output parameter space X is a collection of positive op- 
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erators F x S B(Ti.) such that J2 X F x = 1. A density operator on 7i 
is a positive operator with trace 1. The basic probabilistic interpre- 
tation of these objects is fixed by the prescription that tv(u>T(F x )) is 
the probability to get the measuring result "x" on systems prepared 
according oj, before passing through the channel T. Finally, we call 
a vector maximally entangled, if it is normalized, and its 

reduced density operator is maximally mixed, i.e., a multiple of 1: 

I)*) = {diuiH)- 1 tr(A) . (1) 

Let us set up the equations describing dense coding and teleporta- 
tion in this language. In both cases, the beginning of each transmis- 
sion is to distribute the parts of an entangled state u between sender 
Alice and receiver Bob. Only then Alice is given the message she is 
supposed to send, which is a quantum state in the case of teleporta- 
tion and a classical value in case of dense coding. She codes this in a 
suitable way, and Bob reconstructs the original message by evaluating 
Alice's signal jointly with his entangled subsystem. For dense coding, 
assume that x £ X is the message given to Alice. She encodes it 
by transforming her entangled system by a channel T x , and sending 
the resulting quantum system to Bob, who measures an observable F 
jointly on Alice's particle and his. The probability for getting y as a 
result is then tr(uj(T x ® id)(F y )), where the "<g>id" expresses the fact 
that no transformation is done to Bob's particle while Alice applies 
T x to hers. If everything works correctly, this expression has to be 1 
for x = y, and otherwise (see eq. @). 

Let us take a similar look at teleportation. Here three quantum 
systems are involved: the entangled pair in state u, and the input 
system given to Alice, in state p. Thus the overall initial state is p®uj. 
Alice measures an observable F on the first two factors, obtaining a 
result x sent to Bob. Bob applies a transformation T x to his particle, 
and makes a final measurement of an observable A of his choice. Thus 
the probability for Alice measuring x and for Bob getting a result "yes" 
on A, is tr(p £g> to){F x ® T X (A)). Note that the tensor symbols in this 
equation refer to different splittings of the system (1 <g) 23 and 12 ® 3, 
respectively). Teleportation is successful, if the overall probability for 
getting A, computed by summing over all possibilities x, is the same 
as for an ideal channel, i.e., tr(pA), as in eq. (Q). 

The only relationship between the Hilbert spaces involved, which 
this description requires, is that the input and output spaces of the 
teleportation line are the same, since the whole teleportation process 
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is equivalent to the identity. In some sense the best results (minimal 
dimension for the Hilbert spaces carrying the entangled state, best 
ratio of achieved capacity to capacity used) are obtained in the special 
case, where all Hilbert spaces have the same dimension d, and exactly 
\X\ = d? signals are distinguished. We call this the tight case, and the 
main Theorem refers only to this case. 

Theorem 1 Let 71 be a d-dimensional Hilbert space (d < oo), and X 
a set of d 2 elements. Consider the following types of objects: 

1. Teleportation schemes consisting of 

• a density operator uj on7i®7i 

• a collection of channels T x : B(7i) —* B(7i), x £ X 

• an observable F x , x £ X on Ti ®TL 

such that, for all density operators p on Ti, and A £ B(Ti): 

Y,^(p®u)(F x ®T x (A))=tT P A. (2) 

xex 

2. Dense coding schemes, consisting of the same objects as a 
Teleportation Scheme, but satisfying, instead of ^), the equation 

tt(uj(T x ®id)(F y )) =6 xy . (3) 

3. Bases of maximally entangled vectors, i.e. families of max- 
imally entangled vectors & x G T~i ® Ti, x S X such that 

= S xy (4) 

4- Bases of unitary operators, i.e., collections of unitary oper- 
ators U x £ B(Ti), x £ X such that 

tT(U*U y )=dS xy . (5) 

5. Unitary depolarizers, i.e., collections of unitary operators 
U x £ B{H), x £ X such that for any A £ B(H): 

U* Au x = d tv(A) 1 . (6) 

X 
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Then, given any object of any one of these types, one can construct 
an object of each of the types, using the following equations: 



lo = with $7 maximally entangled. (7) 

F x = | | (8) 

T X (A) = U* X AU X (9) 

® x = (u x ®i)n. (io) 

The logical structure of this result is maybe slightly unusual, so 
we begin by giving some examples how it is used. We can use it, 
for example, as a construction procedure: once we are given a unitary 
basis, we can get from the equations (|7|) to (10) a teleportation scheme 



and a dense coding scheme. Moreover, since we could also start with 
these schemes, ending up with the unitary basis we are assured that 
every teleportation or dense coding scheme is obtained in this way, 
i.e., this construction is exhaustive. In particular, we learn that any 
tight teleportation scheme is necessarily of a very special form: the 
entangled state uj must be pure and maximally entangled, the channels 
T x must be unitarily implemented, and the observable F must be a 
complete von Neumann measurement. 

Another result contained in this Theorem is the amazing equiva- 
lence between (1) and (2): any teleportation scheme works as a dense 
coding scheme and conversely. Alice and Bob merely have to swap 
their equipment to convert one into the other. We must emphasize, 
however, that the tightness condition is absolutely crucial for this 
equivalence. For simplicity, we will discuss this only in the case that 
\X\ = n is not fixed to be d 2 , leaving aside the more difficult question 
what kind of trade-off between resources becomes possible, when u 
lives on H\ <g) Ti.2, with dimensions other than d® d. 

The basic difference between teleportation and dense coding is that 
the parameters d and n have opposite roles: For teleportation d de- 
scribes the size of the signal to be sent, and n describes a resource, so 
the problem becomes more difficult when we increase d and decrease 
n. For dense coding, it is exactly the opposite. Therefore, it is easy to 
show that teleportation (resp. dense coding) schemes exist whenever 
n > d 2 (resp. n < d 2 ). In fact, for teleportation one can take X to 
be a continuum, and replace the sum in the teleportation equation 
by an integral |BD| ], but the dense coding equation would make no 
sense then. The optimality of these dimension inequalities, i.e., that 
no teleportation (resp. dense coding) scheme exists with n < d 2 (resp. 
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n > d 2 ), is also a corollary of Theorem 1. To prove it, suppose we 
had a teleportation scheme with n < d 2 . Then we could add (d 2 — n) 
irrelevant classical signals happening with probability zero (F x = 0), 
and apply the Theorem, which says that all F x must be non-zero after 
all. The same reasoning works for dense coding with the operation of 
throwing in a few unused Hilbert space dimensions. 

Of course, our Theorem is efficient as a construction procedure for 
dense coding and teleportation schemes only to the extent that unitary 
bases can be generated. After giving the proof of the Theorem, we 
will therefore describe the most general construction for such bases 
known to us. 



3 Proof of Theorem 1 
Proof of the Implications "3<^=^4" 

Implicit in the formulation of the Theorem is the claim that the equa- 



tion (10) <& x = (U x (8> not only determines <J> X in terms of U x but 
also, conversely, determines U x in terms of & x . This connection is 
based on a general construction, by which the d 2 matrix elements of 
an operator A : H — » TL are identified with the d 2 components of a 
vector This identification depends on the choice of a maximally 
entangled vector f2. By choosing appropriate orthonormal bases e&, 
k = 1, . . . ,d in the first and second tensor factor, such a vector can 
be written in "Schmidt form" as 

^ = -?= Y] e k <g> e k ■ (11) 
^ d k 

Then a one-to-one correspondence between operators A S BiTL) and 
Vl/ G H <g) H is given by the equation (ek\Aeg) = Vd (ek e||*). We 
will use this in the form 

* = (A® l)tt = (I® A T )n , (12) 

where the transpose operation A t— > A T is defined in the basis e&. 
Then if A and and, similarly, A' and \&' are related in this way, 

(*|(S®I)* , > = 2 tr ( A * BA ') i ( 13 ) 

for arbitrary S G B(7i). Thus ^ is maximally entangled iff this expres- 
sion (for A = A') is equal to <i _1 tr(B), i.e., iff ^4 is unitary. Moreover, 
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setting B = I, the scalar product of vectors ^' is translated to 
d^ 1 tr(A*A') in terms of A, A'. Taking all this together, we get the 
one-to one correspondence between unitary bases and bases of maxi- 
mally entangled vectors, as claimed. Note, however, that this corre- 
spondence depends on the choice of the reference maximally entangled 
vector Q. 



Proof of the Implications "4<^=^5" 

This proof is relatively straightforward, since we are talking about 
only one type of objects, collections of d 2 unitaries U x £ B{TL). It is, 
however, also a crucial step for the entire proof, since it is here that 
the consequences of the tightness condition are seen. We will prove 
this in a form, which is also needed later to establish that the state 
ijj in teleportation and dense coding schemes is necessarily maximally 
entangled. 

The basic observation concerning matching dimensions is the fol- 
lowing. 

Lemma 2 D vectors ,4>d in a, D- dimensional Hilbert space 

form an orthonormal basis if and only if 

D 

£ \<p k )(<p k \ = 1 . (14) 
fc=i 

Of course, this is false when there are more vectors than the di- 
mension of the Hilbert space. Such families of vectors are called "over- 
complete". They exist and are an interesting mathematical structure 
of their own. On the other hand, fewer vectors than the dimension 
can never satisfy eq. (|14|), because the rank (dimension of the range) 
of the operator on the left hand side is at most the number of vectors. 



Proof: It is a well-known fact that eq. (|T4|) holds for any orthonormal 
basis. Conversely, we find from eq. (14) that, for each k, \4>k)(4>k\ < I, 
which is the same as ||</>fc|| 2 < 1- On the other hand, taking the trace 
of (|T|), we get J2k WfikW 2 = tr(I) = D. This is only possible, when 
||(/>fc|| 2 = 1 for all k. Hence the operators \4>k){(t)k\ are hermitian pro- 
jections, and we can invoke the observation that hermitian projections 
PI1P2 with pi +P2 < I are necessarily orthogonal. (For a quick proof, 
sandwich the inequality between factors p±, finding p\ + P1P2P1 < Pi, 
i.e., P1P2P1 = {P2P1)* (P2P1) < 0, and hence p 2 pi = 0). □ 
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We now apply this Lemma to a collection of D = d 2 operators in 
B(Ti), where this space is considered as a Hilbert space with a suitable 
scalar product. 

Proposition 3 Consider d 2 operators K%,... ,K d 2 on a d-dimen- 
sional Hilbert space 7i, and let R > be an invertible operator on 

n. 

Then the following conditions are equivalent 
1. tr(/q/r %) = 5 xy , forx,y = l,...,d 2 
2- E x K*CK X = tr (RC)1 for all C G B(H). 

Proof: Let us define a scalar product {-\-)r on BiTL) by 

(A\B) R = tv(A*R- 1 B) . (15) 

Since R is positive and invertible, this is indeed a scalar product, 
satisfying (A\A) R = only for A = 0. Condition 1 then simply says 
that the K x are an orthonormal basis. By the previous Lemma this is 
equivalent to the completeness relation (|T^) , so all we have to do is to 
show that this relation, adapted to the special scalar product at hand, 
is equivalent to Condition 2 of the present Lemma. The completeness 
relation is that, for any A, B G B(Tt), 

(A\B) R = ^2(A\K X ) R (K X \B) R . (*) 

X 

It suffices to evaluate this on rank one operators A, B G B(H), since 
these span the whole space. We take A = \4>i)(4>2\ and B = IV'iXV^I- 
Then the left hand side of equation (*) becomes 

(<£l|irVl)W*> , (*LHS) 
whereas the right hand side is 

Yj^AR^K^^KRT 1 ^) = (V2|M0 2 ) , (*RHS) 

X 

with 

m = Y,KR~ 1 \^i\R' 1 k x = J2KCK X 

X X 

where we have interchanged the two factors in each term, and intro- 
duced the abbreviation C. Since (*LHS)=(*RHS) for every tp2, <t>2-, we 
find M = 1. The factor is readily identified as 
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(<^i|i? _1 V>i) = tr(RC). Since operators of the form C span B(Tt), 
the completeness relation thus becomes equivalent to ^2 X K*CK X = 
tr (RC)TL for all C, which completes the proof. □ 

The special case of this Proposition, where each K x is unitary 
and R = 4 1, is exactly the relationship between items 4 and 5 of 
Theorem 1. However, there is another consequence needed later on: 

Corollary 4 Let 17%,. . . , U^p. £ B(TC) be unitaries in a d- dimensional 
Hilbert space 7i, and p a density operator such that ti(U*pU y ) = 5 xy . 
Then p = d TL. 

Proof: Since the U x are an orthonormal set whose cardinality is the 
dimension, there can be no null vectors of this scalar product, i.e., 
tr(A*pA) = implies A = 0. Hence p is invertible, and we can apply 
the previous Proposition with R = p^ 1 , finding that ^}2 X U*AU X = 
tr(p -1 A)I. The trace of this equation is d 2 ti(A) = dtr(p' 1 A). This 
holds for all A, i.e., p' 1 = dl. □ 

Proof of "(3 or 4)=^(1 and 2)" 

Suppose now we are given either a basis of unitary operators or of 
maximally entangled vectors. Then we can choose a maximally en- 
tangled vector f2 and use equation (|l0|) as in in the proof of "3<^>4" 
to define the other kind of basis. Equations (Q) and (|9|) then become 
explicit definitions of the observable F x and the transformations T x , 
respectively, so all the objects needed for a teleportation or dense cod- 
ing scheme are defined, and we only need to verify that equations @ 
and (H|) are indeed satisfied. 

In the teleportation equation an expectation value is generated 
between a state on the first and an observable on the third factor of a 
triple tensor product. This is a consequence of a similar "teleportation 
equation" on the level of vectors, which we now state. For later use 
we prove a certain converse at the same time. 

Lemma 5 Let E C d (g> C d be the maximally entangled vector Q, = 
d^ 1 / 2 J2k e k ® &k, where e^, k = 1, . . . , d is the standard basis of C d . 
Let M G B(C d ), and /i6C. Then the equation 

((f) <g> <g> m i)n (8> v) = /i(</#) 
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holds for all (j),ip G C d , if and only if M = dpi. 
Proof: Inserting the sum denning £1 we get 
(0® n|(I® M® 1)Q ® V>) = 

= i^(0(g>e fc ® e fc |(l® M® I)e^®e^®V) 

= ~£(0M(e*|Me^)(e fe |</>) = i(0|M T ^} , 

which is equal to p(<^\tp) for all (f),ip iff M T = dpi. □ 



Consider now the term with index x G X in the teleportation 
equation (g), with and T x defined via equations (||), (|9|), and (10). 



Without loss of generality we set p = |0i)(02|, A = IV'iXV^I- Then 

term x = (0 2 ® ft|$ x ® U*^x){^ x ® f/*V>2|0l ® ^) • 

The first scalar product can be rewritten by substituting $ x from 
equation (10), using equation (|l2|): 



(0 2 ® n\$ x ® C/*^i) = (0 2 ® ((1® ^)^)|((£4® J)^)®^i) 
= (02® ((C/J® I)0)|((I® t/J)0)®Vi) 
= <(l®?/J®I)^2®n|(l®Z7j®I)f2®^i) 

= (0 2 ®ft|Sl®Vi} = -;(02|^i) , 

where at the last equation we used Lemma [5] with p = 1/d. To- 
gether with a similar computation for the second scalar product, we 
get terra,; = d _2 (02|V'i)(^2|0i) = d~ 2 tr(pA), and equation @ follows 
by summing over d 2 equal terms. 

Similarly, for the dense coding equation (||) we get 

tT(u(r x ®id)(F v )) = (n\(u* ® i)$ v )($ v \(u x ® i)n) , 

i.e., the absolute square of the scalar product 



1 



(n\(u; i)$ y ) = (n\(u*u y ® i)n) = -ti(u* x u y ) = s x 



Jxy ; 



where we have used, in turn equation (|10[), the maximal entangled- 
ness of f2 (see equation (|l|)), and the orthogonality of the U x . This 
completes the proof of the dense coding property. 
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Proof of the Implications "2 =4> Rest" 



Let us now assume that a dense coding scheme is given. We have to 
conclude that it is of the special form given in equations (|7|...|l0|). 

Note first that if u = J2a (A Q > 0) is a mixture of states 

satisfying the teleportation equation, then every tu a also satisfies it. 
Hence the assumption is also satisfied for each pure component u a , 
and we can first analyze the problem assuming uj to be pure. In order 
to show that uj indeed is pure, we only have to verify that the given 
F, T are consistent only with one pure state. So for the moment we 
will assume that uj = is pure. 

The next step is a simple general observation on the coding of clas- 
sical information on quantum channels, which we isolate in a Lemma. 

Lemma 6 Let K be a D- dimensional Hilbert space, and cr x ,F x £ 
B(fC), for x £ X , a set with D elements. Suppose that each a x 
is a density operator, F is an observable, and tr(a x F y ) = 5 xy , for 
x,y = l,... ,D. 

Then there is an orthonormal basis & x £ TC such that 

a x = F x = \<S> X )(<S> X \ . 

Proof: Let x be one of the normalized eigenvectors of o~ x with non- 
zero eigenvalue. Then since F x < 1, and ($> X \F X $ X ) = 1, & x must also 
be an eigenvector of F x with eigenvalue 1. Similarly, for any y ^ x the 
F x > 0, and the normalization J2x F X = ~R. forces F y <fr x = 0. Hence the 
$ x are orthonormal, and since their number is the dimension of the 
space, they must be a basis. Consequently we have jointly diagonalized 
the F x and the a x , with eigenvalues either or 1. □ 

We apply this Lemma with D = d 2 and a x the state after applica- 
tion of T x to the first factor, i.e., ti(o~ x A) = tr(a;(T2; (g) id)(A)). This 
proves equation (Js|) , although it remains to be seen that each & x is 
maximally entangled. 

Since the a x form a maximal set of pure states, there cannot be a 
non-zero projection P such that, for all x £ X, 

= tr(a x (l® P)) = tr(uj{T x (g>id){l(g> P)) 

= tr(o;(i® p)) = p)n) . 

Hence must have full Schmidt rank. We will need the consequence 
that the equation (A 1)0 = (A' implies A = A'. 
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Let T X (A) = J2 a Kx a AKx,a be the Kraus decomposition of T x . 
Then the teleportation equation is 

a 

Therefore, {{K xa (8> 1)01$^) = for all y ^ x, and for every x there 
must be constants c a such that 

(K x , a ® 1)0 = c a $ x . (16) 

Since O has full Schmidt rank, this implies that all K XtQ . are propor- 
tional to each other, i.e., that T x can be written with a single Kraus 
summand. Of course, the corresponding K x = U x must be unitary, 
and since both sides are normalized, equation ( |l6| ) & x = (U x (g> 1)0, 
possibly after fixing suitable phase factors (which influence neither T x 
nor F x ). 

The orthonormality of the Q x translates into tr(pU*U y ) = 8 xy , 
where p is the reduced density operator of oj. But then Corollary |] 
shows that p must be a multiple of the identity, i.e., O and each Q x is 
maximally entangled. 

Finally, we have to complete the argument for the purity of uj by 
showing that only one pure state is consistent with the other data 
T,F, encoded in U x . But this is obvious from the explicit expression 
= (U* ® l)$ x . 

Proof of the Implications "1 =^> Rest" 

Let us now assume that a teleportation scheme is given. We have to 
conclude that it is of the special form given in equations (|7t..|l0|). 

The crucial input for this proof is the principle that in quantum 
mechanics there is no measurement without perturbation. It enters 
in the following form, a corollary of the so-called Radon-Nikodym 
Theorem for completely positive maps. We state it here as a Lemma. 

Lemma 7 Let 7i be a finite dimensional Hilbert space, and let T a : 
B(7i) — > B(7i) be completely positive maps such that ^ a T a = id. 
Then there are positive numbers t a such that T a = t a \&. 

Proof: For readers less familiar with dilation theory of cp-maps we 
include a quick proof based on the Kraus decomposition T{A) = 
J2p KpAKp, which exists for every completely positive map. Note 
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that by decomposing each T a in Kraus form, we get a finer decompo- 
sition of id, so we may as well prove the Lemma for the case that each 
T a is of the form T a (A) = K*AK a . With A = |^)(^|, 

a. 

Hence K*tp = X(ip)ip, with a factor £ C. But then every vector 
tp is an eigenvector of the linear operator K^, which is only possible, 
if K* is a multiple of the identity. □ 



A collection of completely positive maps adding up to a normal- 
ized one should be understood as an "instrument" in the terminol- 



ogy of Davies [Da], i.e., a device which produces classical measure- 
ment results "A;", such that the probability for obtaining this result 
and a response to a subsequent measurement F on an input state p 
is tr(pTfc(F))). The channel J2k^k then describes the overall state 
change, when the measuring results are ignored. In this language the 
hypothesis of the Lemma says that there is no overall state change 
through the device, i.e., "no perturbation" of the system. The con- 
clusion is that in that case the output probabilities are tk, and inde- 
pendent of the input state, i.e., no information about the system is 
obtained. 

As a first application, we conclude exactly as in the previous sub- 
section that each convex component of the state u again satisfies 
the teleportation equation. Hence we can once more assume that 
uj = |fi)(fl| is a pure state. The argument that 0, is then uniquely 
determined by the other data, and hence that u is pure is the same 
as in the dense coding case. 

Clearly, this kind of argument is also useful for decompositions 
of T x or F x into sums of (completely) positive terms. To do this 
systematically, fix a maximally entangled unit vector 3, so that vectors 
m.7i®7i become expressed as $ = (A(g)I)3 for a uniquely determined 



operator A (see equation (12)). In particular, we can write Q = (PF(g> 
1)3, and the Kraus decomposition and spectral decomposition of each 
F x in the form 

T X (A) = J2K* AK *,<x ( 17 ) 

a 

F x = E(^® I) 1^1(^0 1)*. (18) 
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Inserting this into the teleportation equation (Q) we find a sum over 
x, a, P, in which each term represents a completely positive operator, 
and which sum up to the identity. Hence by Lemma |7|, each term has 
to be multiple of the identity, Jl x a /jid, say. This can be written in 
terms of scalar products, if we take p = |0i}(02| anci A = iV'iXV^h 

V>x,a,P tr(p A) = Jl x> a,l3 (<^2,'0l}(^2,0l) 

= (<f>2 <8 Q, (Ar l/3 (8) I (8) i^, a )3 <g> Vi) 

(3 (8) ^2, 8) 1(8) if a; ,a)01 ® fi) . 

Note that the two scalar products on the right hand side are com- 
plex conjugates of each other apart from a swapping of the arguments 
(<fe) ^l) an d (^2)01)) which exactly matches the variable pairing on 
the left hand side. Since the equation is to hold for arbitrary vectors 
01 , 02 ) ipi ) "02 ) we can hold one pair fixed and find that 

<0 2 <8 n, (A XiP ® I <g> if* JE ® Vi) = Mz,a,/3 (02, ^1) , (19) 
where p x ,a,p is a factor determined in terms of ju, and the scalar prod- 



ucts involving (?/>2,0i). With 17 = (W (8 1)3, and equation (12) we 
get 

<0 2 8> n, (A Xi p ® i <8 isr* a )H ® fa) 

= (02 ® (I <8 1^,0)3, l)((^,/3 8) 1)3 (8 Vi)) 

= (02 (8 (ifJ Q (8 1)3, (I® W* (8 I) ((I (8 A^)H (8 ^1)) 
= (0 2 (8 3,(I(8X^PF*A^ / g(8l)3(8'0i) 



(02,^ 



where we have used the notation K = (K*) T for the matrix element- 
wise complex conjugation in the Schmidt basis belonging to the max- 
imally entangled state H. Since the above equation holds for all 02 
and ipi ) Lemma [| implies that 



K X)Cl W*A^ = d nv^p 1 , (20) 

for all x, a, (5. 

Let us say that a label x £ X contributes to teleportation, if the 
corresponding term in the teleportation equation does not vanish for 
all p and A. This is equivalent to saying that for some a, (3 the factor 
Ma; a/3 is non-zero. For such triples (x, a, (3) all three operators on the 
left hand side of equation (|(J) have to be invertible. 
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Now since there has to be at least one contributing label, W has 
to be non-singular, which means that £1 has full Schmidt rank. Equiv- 
alently, the reduced density operator u\ for the first factor has no zero 
eigenvalues. From this we conclude that the non-contributing labels 
are precisely those for which F x = 0. Indeed, we may set A = p = 1, 
and use the normalization of T x to find 

= tr((I®w)(F«® I)) =tr((I®wi)i^) 

Since F x > 0, and I (g) u\ has only strictly positive eigenvalues, this 
implies F x = 0. 

Now let x be a contributing index, and choose some triple (x, a, (3) 
with (i x ,a,p 7^ 0- If we now a t equation (|20|) for triples (x, c/,/3) 
with arbitrary a', we get if X)a ,' = (p> x ,a',p/ ' Hx,a,p)Kx,a, i-e., all Kraus 
operators of T x are proportional, and hence T x can be written with a 
single Kraus summand, T X (A) = U*AU X , with a unitary f^. 

Similarly, we find that all A x ^/ are proportional, which means that 
F x = I^X^xl with ^ = (A s ®'l)H. 

We can now apply Lemma |2| to these vectors & x , setting & x = 
for non-contributing labels. The conclusion is that the & x are an 
orthonormal basis. In particular, all indices do contribute after all. 

Equation (^C|) and the unitarity of U x allow us to express A x in 
terms of U x : 

A x = dfi x U x W~ l (21) 
Orthonormality of the $ x becomes 

5 xy = i tr(^*^) = tr(^ W'iW 1 )* U*) . (22) 

For x = y we find that |^a;| 2 is independent of x, hence the operators 
(jXx~/\n x \)XJ* are unitary, and satisfy the hypothesis of Corollary |] with 
p a positive multiple of W (W )*. Hence this operator is a multiple 
of the identity, W is unitary up to a factor, and O = (W (8> I)S is 
maximally entangled. Moreover, we see from equation ([22] ) and the 
U x form a unitary basis. 

Since H was an arbitrary maximally entangled vector, we may just 



as well take 5 = so equation (21) holds with W = I. Hence, 
3" is = c (fj,; (8> I)f2, where c is a factor which has to be of modulus 1, 
because f2 and $ x are normalized, and U x is unitary, and which can 
be chosen to be 1 by adjusting the phase of & x . This completes the 
proof. 
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4 Constructing bases of unitaries 



It is not a priori clear that bases of unitary operators should exist 
in any dimension. Indeed, the system equation (|5|) of equations is 
formally overdetermined, according to the following rough dimension 
count. The variables in this system are the unitaries U x , each of 
which we can take in the (d 2 — l)-dimensional manifold SUd, i.e., with 
det(U x ) = 1, by fixing a phase factor. Since the transformations U x \— > 
ViU x V2, for arbitrary Vi, V2 £ SUd leave the set of solutions invariant, 
we may fix U\ = I, and take U% diagonal without loss of generality. 
This reduces the number of variables to (d — 1) + (d 2 — 2)(d 2 — 1). 
On the other hand, orthogonality introduces one complex constraint 
for every pair x 7^ y. None of these is trivially satisfied due to the 
special choices we made, so we have to take d 2 (d 2 — 1) constraints into 
account. This leaves, formally, 

#variables-#equations = — [d — l)(2d + 1) < . 

Of course, we know that this count is somehow too crude, because, 
after all, many inequivalent unitary bases are constructed below. But 
it is not so easy to spot the dependences among the constraints. Note 
also that the dimension count is essentially the same for bases orthog- 
onal with respect to a weight p 7^ d~ 1, but in that case Corollary |] 
shows that there is no solution at all. 

In order to describe the best known construction for unitary bases 



[VW], let us introduce some terminology. We say that a (single) uni- 
tary matrix is of shift and multiply type, if it is the product of a 
permutation operator and a diagonal unitary. In other words, every 
row or column contains {d— 1) zero entries, and one entry of modulus 
1. The bases we will construct not only have the property that each 
element is of this type, but also that the d 2 values for x can be split 
into d options for "shift" and d options for "multiply". 

Definiton 8 A shift and multiply basis of unitary matrices in C d 
is a collection of d 2 unitary operators Uj, i,j G Id = {1, • • • ,d}, 
satisfying the orthogonality relation ti{U*jUki) = d 5ikd~j£, and acting 
on the basis vectors \k) as 

U ij \k)=Hf k \X(j,k)) , (23) 
where the Hj k are complex numbers, and A : Id X Id — > Id- 
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Proposition 9 The parameters and A : Lj x Lj — > Id define a shift 
and multiply basis of unitary matrices if and only if the following two 
conditions are satisfied 

1. Each ffl is a Hadamard matrix, i.e. \Hf k \ = 1 for all i,k, 
and W{W)* = d I. 

2. X is a Latin square, i.e., the maps k i-> \(k,£) and k \(£, k) 
are injective for every I. 

Proof: For Uij to be unitary, it is necessary and sufficient that the Hf k 
are phases, and that k i— > X(j, k) is injective (hence bijective) for every 
j. For the orthogonality we have to evaluate 

tr(t/*.c/ iY ) = YMk H k(KmWM . 
k 

We consider first the case j = j'. Then the scalar products in the sum 
are all equal to 1, and equating this expression to Stf we find that H 3 
is Hadamard. 

Now let j 7^ j', and consider the "coincidence set" C = {k | 
X(j,k) = \(j',k)}. Then orthogonality requires, for every i,i', that 

= E = E < Xo(fc) (^Jw = (Hi'xcH^U , (24) 

fcec fe=i 

where xc(^) = 1 f° r k £ C, and zero otherwise, and in the last line xc 
denotes the projection xc\k) = Xc(k)\k). But since H 3 and H 3 are 
Hadamard, and in particular invertible, this implies \c = 0. Hence C 
is empty, and the second injectivity of A is proved. □ 

In order to construct unitary bases of this form, we must now 
construct Hadamard matrices and Latin squares of the appropriate 
dimension. For both of these tasks there is a rich literature, and 
below we will give a brief summary on what is known for each. 

It is useful to note that each of the structures 'unitary bases', 
'Hadamard matrices', and Latin squares' has a natural notion of 
equivalence, and to some extent these equivalences are related. We 
call two unitary bases U,U' equivalent, if U' x = V\U x iV2, for some 
unitaries V\,V2, and a re-labelling x i— ► x' . Hadamard matrices are 
called equivalent, if one is obtained from the other by permuting rows 
or columns, or multiplying rows or columns with phases. Finally, a 
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Latin square A : Id x Id — > Id is equivalent to any other obtained by 
applying a permutation on each of the three copies of Id involved. In 
each case there are also discrete transformations, such as transposi- 
tion or complex conjugation (where applicable). It should be noted 
that replacing each W by an equivalent one, typically only leads to 
an equivalent unitary basis, if the equivalence operation is the same 
for each j. With j-dependent equivalence transformations it is pos- 
sible to construct inequivalent unitary bases in d = 3, although in 
this dimension there is only one Hadamard matrix and only one Latin 
square - up to equivalence. Of course, in d = 2 all three structures, 



including the unitary bases are unique up to equivalence | V"W |. The 
unique unitary basis is then given by the three Pauli matrices and the 
identity and, of course generates via the Theorem 1, the usual two 
qubit examples of teleportation and dense coding. 

For each of the three structures we furthermore have an obvious 
notion of tensor product, allowing the construction of a unitary basis 
(resp. a Hadamard matrix, or Latin square) in dimension d = d\d2, if 
counterparts in dimension d\ and dimension e?2 are given. 

In order to show that unitary bases exist in any dimension it is 
easiest to use group theory based constructions: the Latin square 
can be taken as the multiplication table of any group of order d, for 
example the cyclic group. The Hadamard matrix can be taken as the 
matrix implementing the Fourier transform on an abelian group of 
order d, the standard example being given once again by the cyclic 
group of order d. Thus H^i = exp(^p/c£), where k and I are taken 
modulo d. If we combine these data into a unitary basis we get an 
instance of what we propose to call a unitary basis of group type ("nice 
error basis" in JKn| | ) . These are orthonormal unitary bases with the 
additional property that the operator product of any two elements is 
a third, up to a phase. That is to say, the index set X is a group, and 

U x U y = n(x,y)U xy , (25) 

with \fi(x,y)\ = 1. In the special case of an abelian group X this 
is a discrete version of Weyl systems of unitary operators, named af- 
ter their continuous variable counterpart, well-known from quantum 
optics and non-relativistic "phase space" quantum mechanics. 

Latin squares are not completely classified, nor does there seem 
to be a realistic hope to do so. A standard work on the subject is 
[DK], a useful net resource is Rj]. Counts of squares are usually done 



for "normalized squares", in which the first row and column are in 
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natural order, thus eliminating some trivial freedom. In d = 5 Euler 
counted 56 of these, but only 2 are inequivalent, because the symbols 
themselves can also be permuted. Counts of normalized squares have 
now gone all the way up to d = 10, but are no longer done by hand 
(there are roughly 7.5 x 10 24 [MR]). It is also clear from these numbers 
that group based constructions exhaust only a tiny fraction of the 
possible unitary bases. 

Hadamard matrices are also a standard subject in coding theory. 
However, usually only the real case (orthogonal matrices with entries 
±1) is considered. It is easy to see that real Hadamard matrices exist 
only in dimension two and multiples of four. Again, the possibilities for 
such designs by far exceed the group based possibilities (the characters 
of an abelian group are real only if d = 2 n ). A standard reference is 



For complex Hadamard matrices the Fourier matrices show that 
there is no constraint on dimension. The uniqueness in d = 3 is easy 
to get. The general form in d = 4 is, up to equivalence 



/ 1 
1 
1 

V 1 



1 



1 \ 

-1 -1 

u —u 
-u u 



(26) 



where u is an arbitrary phase. For u = 1 this is equivalent to the 
Fourier matrix of "Klein's Four Group", the product of two copies of 
the two-element group, and for u = i it is equivalent to the Fourier ma- 
trix of the cyclic group. The possibility of embedding the cyclic group 
Fourier matrix into a higher dimensional manifold can be generalized 
to arbitrary composite numbers d = pq: whenever Vke is a matrix of 
phases satisfying the periodicity conditions Vk,i = Vk+ Pt £ = Vk/+ q , we 
get a Hadamard matrix as 



'2ni 

Hke = V k e exp ( — hi 



(27) 



One might conjecture from this that for prime orders d the Hada- 



mard matrix is unique. This problem was discussed by Haagerup Ha ] 
on the basis of a completely different motivation (theory of von Neu- 
mann algebras). There it is shown that d = 5 there is uniqueness, but 
for d = 7 there are at least 5 solutions. For some primes, uncountably 
many inequivalent Hadamard matrices are known. 
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